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A series of three-dimensional numerical simulations is used to study the intrinsic stability of high-speed 
turbulent flames. Calculations model the interaction of a fully-resolved premixed flame with a highly subsonic, 
statistically steady, homogeneous, isotropic turbulence. The computational domain is unconfined to prevent 
the onset of thermoacoustic instabilities. We consider a wide range of turbulent intensities and system sizes, 
corresponding to the Damkdhler numbers Da = 0.1 —6.0. These calculations show that turbulent flames in the 
regimes considered are intrinsically unstable. In particular, we find three effects. 1) Turbulent flame speed, 
St, develops pulsations with the observed peak-to-peak amplitude > 10 and a characteristic 

time scale close to a large-scale eddy turnover time. Such variability is caused by the interplay between 
turbulence, which continuously creates the flame surface, and highly intermittent flame collisions, which 
consume the flame surface. 2) Unstable burning results in the periodic pressure build-up and the formation 
of pressure waves or shocks, when St approaches or exceeds the speed of a Chapman-Jouguet deflagration. 
3) Coupling of pressure gradients formed during pulsations with density gradients across the flame leads to 
the anisotropic amplification of turbulence inside the flame volume and flame acceleration. Such process, 
which is driven by the baroclinic term in the vorticity transport equation, is a reacting-flow analog of the 
mechanism underlying the Richtmyer-Meshkov instability. With the increase in turbulent intensity, the limit- 
cycle instability discussed here transitions to the regime described in our previous work, in which the growth 
of St becomes unbounded and produces a detonation. 


I. INTRODUCTION 

Exothermic reaction fronts play a fundamental role in 
systems ranging from the propulsion and energy con¬ 
version applications on Earth, e.g., automotive and air¬ 
craft engines or power-generation turbines^, to the as- 
trophysical thermonuclear explosions of the white dwarf 
stars (Type la supernovae), which are the key sites for 
the galactic nucleosynthesis and production of heavy ele¬ 
ments in the Universe.^ The dynamics of burning in such 
systems is critically controlled by a rich variety of un¬ 
stable phenomena exhibited by flames. In the context of 
premixed flames, which are the focus of this study, most 
notable examples include hydrodynamic (e.g., Landau- 
Darrieus), various thermodiffusive (e.g., cellular), and 
body-force (e.g., Rayleigh-Taylor) instabilities in laminar 
flames.^ These are intrinsic flame instabilities, as their 
onset is not dependent on the external factors, such as 
the details of the upstream flow or the overall geometry 
of the combustion system. 

In contrast, instabilities of turbulent flames are pri¬ 
marily considered in the context of burning in confined 
environments with walls or obstacles. They are viewed to 
result from the resonant coupling between the exother¬ 
mic process and the acoustic field, which it generates in 
the interior of a combustor.Thus, such thermoacous¬ 
tic instabilities cannot be considered in isolation from the 
host system. 

Thermoacoustic instabilities of turbulent flames have 
been studied extensively using experimental, theoretical, 
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and numerical techniques.At the same time, a more 
basic question of the intrinsic stability of unconfined, pre¬ 
mixed turbulent flames remains largely unexplored. It is 
unclear whether a flame propagating into a statistically 
steady, homogeneous, isotropic upstream turbulence in 
the absence of the surrounding combustor walls can ex¬ 
hibit significant variations in its key dynamical charac¬ 
teristics. Of course, even in such an idealized situation, 
turbulent flame speed cannot be exactly constant due to 
the inherent variability of the turbulent flow field. The 
question, however, is: In a statistically steady turbulent 
flow, can the state of a flame be meaningfully described 
by a single value of the flame speed as is typically as¬ 
sumed in modern turbulent combustion models^’®’®, or 
whether the burning rate can exhibit variations with the 
magnitude comparable to its mean value or even larger? 

In realistic combustion systems, the large-scale flow 
is rarely homogeneous, isotropic, or statistically steady. 
Modern numerical models of such systems, e.g., a ram¬ 
jet combustor^*^ or the white dwarf interior in a Type 
la supernova explosion^ typically follow the large-eddy 
simulation (LES) approach and resolve the unsteady 
large-scale flow explicitly.^’® At the same time, they cru¬ 
cially rely on turbulent combustion models to capture 
the dynamics on sub-filter scales, where the flow is lo¬ 
cally isotropic, homogeneous, and statistically steady rel¬ 
ative to the characteristic dynamical time of large scales. 
While the intrinsic instability of a turbulent flame, if 
present, would affect the flow on all scales, it is the small 
sub-filter scales, where the question posed above is par¬ 
ticularly relevant. 

Here we study the intrinsic stability of turbulent flames 
in the regimes characterized by the Damkohler number 
Da = Tf^d/TR ~ 1 and Karlovitz number Ka = tr/tr > 
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1 , where Ted, tji, tk are, respectively, the characteris¬ 
tic large-scale turbulent, reaction, and Kolmogorov time 
scales (“thin reaction zones regime”®). 

These regimes were chosen for two reasons. First, ear¬ 
lier calculations by Poludnenko and Oran^^’^® of the in¬ 
teraction of a premixed flame with a statistically steady, 
homogeneous, isotropic turbulence showed variability of 
the turbulent burning speed, St, with the maximum 
peak-to-peak amplitude « 3. Those stud¬ 

ies considered turbulence with the integral velocity Ui = 
18.5S'l and scale I = 1.9i5l, where Sl and 6l are, re¬ 
spectively, the speed and the thermal width of a laminar 
flame. In such high-speed turbulence, the characteristic 
turbulent time scale is much shorter than the reaction 
time scale, i.e., Da 1. At the same time, at lower tur¬ 
bulent intensities or in larger systems, i.e., when Da ^ 1, 
Ted and Tfl become comparable. While turbulence acts 
as a perturbing force increasing the flame surface and 
the overall burning rate, flame propagation acts as a 
restoring force consuming the flame surface and decreas¬ 
ing the flame speed. Thus, in the regimes characterized 
by Ted ~ a resonant state may develop between the 
turbulence and the flame, which could result in a larger 
variability of the burning speed than observed by Polud¬ 
nenko and Oran^^ds 

Second, previously Poludnenko et al}'^ showed that 
flame propagation is not possible in sufficiently fast tur¬ 
bulence with Da < 1 and Ka^ 1. Above a certain sub¬ 
sonic turbulent intensity and system size, an unconfined 
flame spontaneously develops a strong shock or a deto¬ 
nation. Thus, it is important to understand the flame 
dynamics in the lower-speed regimes below this critical 
stability threshold, in which a turbulent flame can still 
exist but its properties may be significantly affected by 
the compressibility effects. 

Variability of the burning speed of a premixed flame in¬ 
teracting with a homogeneous, isotropic upstream turbu¬ 
lence can also be seen in the direct numerical simulations 
(DNS) by Nishiki et 17_ contrast 

to Poludnenko and Oran^^ds^ those studies considered 
much lower intensity turbulence with Da « 17 — 18 in 
Nishiki et and Da « 1.5 in Bell et , with 

Ka ^ 1 in both cases. Nevertheless, the maximum 
observed peak-to-peak amplitude < 2.0 is 

comparable to that found by Poludnenko and Oran^^ds 
Those calculations, however, had two important limita¬ 
tions. Nishiki et a 7 ^®d 6 followed the flame evolution only 
for a relatively short time after the flame became fully de¬ 
veloped, namely l.bTed, compared to in Polud¬ 

nenko and Oran^^ds Calculations of Bell et , on 
the other hand, lasted for ^Ted- However, they were per¬ 
formed in the two-dimensional geometry and, thus, did 
not have a realistic turbulent flow field. Nevertheless, 
results of those studies also suggest that flame dynam¬ 
ics can be quite unsteady even under the most idealized 
conditions. 

Other prior three-dimensional (3D) DNS studying sim¬ 
ilar turbulent combustion regimes {Da < 1, Ka > 1) gen¬ 


erally do not allow the question of the intrinsic stability 
of turbulent flames to be addressed. They either con¬ 
sider time-varying turbulence, e.g., decaying^® or shear- 
driven^*^, study flame configurations that inherently re¬ 
sult in a time-varying flame evolution, e.g., spherically 
expanding turbulent flames^ or have domain sizes with 
unrealistic turbulent integral scales that are smaller than 
a characteristic laminar flame width^^ (also see Ref. [23]). 

In this paper, we present a systematic study of the 
intrinsic stability of premixed turbulent flames using 
a series of fully-resolved 3D calculations. In order to 
eliminate any potential external sources of instability, 
we consider a globally unstrained, statistically planar 
flame interacting with a statistically steady, homoge¬ 
neous, isotropic turbulence in the absence of external 
walls or boundaries. Calculations described here model 
a broad range of turbulent conditions and reacting mix¬ 
ture properties. Special emphasis is made on following 
the flame evolution over a large number of characteris¬ 
tic eddy turnover times to demonstrate the long-term 
behavior. Results of these calculations show that, in a 
certain range of regimes, premixed turbulent flames are 
inherently unstable. 

II. NUMERICAL CALCULATIONS 
A. Model and method 

Turbulence-flame interactions are modeled using com¬ 
pressible reacting-flow equations with molecular trans¬ 
port processes and chemical energy release 
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Here p is the mass density, U is the velocity, E is the 
energy density, P is the pressure, Y is the mass fraction 
of the reactants, q is the chemical energy release, and Y 
is the reaction source term. The coefficients of thermal 
conduction, V, and molecular diffusion, D, are 

D = DqT^, 1C = KoCpT^. (5) 

H is the viscous stress tensor 

n = M(vu-f (vu)T)-| mv-u)i, (6) 

where I is the unit tensor, and the coefficient of shear 
viscosity is 


p = MoT". 


(7) 
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In eqs. (5) and (7), Dq, kq, /io, and n are constants, 
and Cp = 7 i?/M (7 — 1) is the specific heat at constant 
pressure. 

The equation of state is that of an ideal gas. Chem¬ 
ical reactions are modeled using first-order, single-step 
Arrhenius kinetics 

dY ■ f Q \ 

— sy = -pyBexp(-—), (8) 

where B is the pre-exponential factor and Q is the acti¬ 
vation energy. 

The above physical model is supplemented with a 
mechanism of turbulence generation in the domain. An 
approach often used in other studies of turbulence-flame 

interaction^®d7,24 jj^volves injecting at the upstream do¬ 
main boundary a predefined turbulent velocity field, 
which is subsequently advected downstream toward the 
flame. The advantage of this approach is that no external 
turbulence forcing is required in the domain. The draw¬ 
back, however, is that as turbulence is being advected 
through the domain, it decays and, thus, in general it 
is time-varying and inhomogeneous. The characteristic 
time scale of turbulence decay is the large-scale eddy 
turnover time, Ted- Therefore, in situations character¬ 
ized by Da > 1, or Ted ^ tr, the rate of this decay is suf¬ 
ficiently slow compared to the reaction time scale and, 
thus, turbulence can be viewed as statistically steady 
on relevant dynamical times. In contrast, this approach 
is not applicable at high turbulent intensities^^, i.e., at 
Da < 1, which are the subject of this study. In these 
regimes, turbulence would decay on time scales shorter 
than the laminar flame self-crossing time, and often be¬ 
fore it can be advected from the upstream boundary to 
the flame. Instead, in high-speed flows, steady-state tur¬ 
bulence requires active driving inside the computational 
domain. 

In the calculations discussed here, such driving, repre¬ 
sented by the forcing term T in eq. (2), is implemented 
using a spectral method.In this approach, velocity 
perturbations (5U(k) are initialized in the Fourier space. 
Each component <5t/i(k) is an independent realization 
of a Gaussian random field superimposed with the de¬ 
sired energy injection spectrum of arbitrary complexity. 
An inverse Fourier transform of 5U(k) gives 5U(x), the 
velocity perturbation field in the physical space. It is 
ensured that: (1) <5U(x) is divergence-free, i.e., turbu¬ 
lence driving does not artificially induce compressions or 
rarefactions; (2) (5U(x) provides a constant rate, e, of 
kinetic-energy injection; and (3) the total momentum in 
the perturbation field is zero so that no net momentum 
is added to the domain. Resulting velocity perturbations 
are added to the flow field U(x) at every time step, and 
the overall perturbation pattern is periodically regener¬ 
ated. Detailed description of the method can be found in 
Poludnenko and Oran^^, and the analysis of the result¬ 
ing turbulence, both reacting and non-reacting, includ¬ 
ing comparison with prior experimental and DNS results, 
was presented in Refs. [27, 28]. 


TABLE I. Input model parameters and resulting computed 
laminar flame properties 


Parameter 

Value 

Cases'^ 

To 

Input 

293 K 


Po 

1.01 X 10® erg/cm^ 


PO 

8.73 X 10“'‘ g/cm^ 


7 

1.17 


M 

21 g/mol 


B 

6.85 X 10^^ cm®/(g s) 



1.71 X 10^^ cmV(g s) 

SlOs, S16s 


1.45 X 10^^ cmV(g s) 

S16sa 

Q 

46.37 RTo 



27.57 RTo 

S16sa 

q 

43.28 RTo / M 



18.20 RTo / M 

S16sa 

Ko = Do 

2.9 X 10“® g/(s cm K") 



7.25 X 10“® g/(s cm K") 

SlOs, S16s 


1.06 X 10“® g/(s cm K") 

S16sa 

Po 

3.1 X 10“® g/(s cm K") 

S15v 

n 

0.7 


Tpfi 

Output 

2135 K 



1068 K 

S16sa 

ppfi 

1.2 X 10“'‘ g/cm® 



2.4 X IQ-'^ g/cm® 

S16sa 

Sl 

0.032 cm 


Sl 

302 cm/s 



75.5 cm/s 

SlOs, S16s, S16sa 


^ Multiple values are shown for parameters that differ between 
calculations. The first line for such parameters gives values 
corresponding to the base reaction-diffusion modeP^’^^, 
subsequent lines give values specific to the calculations 
indicated in the rightmost column. 


Table I lists values of various parameters of the physi¬ 
cal model described by eqs. (1) - (8). The default values 
of the parameters are calibrated to reproduce the cor¬ 
rect laminar flame properties of a stoichiometric H 2 -air 
mixture under Le = 1 conditions. Details of this base 
reaction-diffusion model are given in Refs. [12, 13, 29]. 
Some model parameters were also varied to study the ef¬ 
fect of the laminar flame Mach number and the density 
ratio across the flame on the flame stability (see § IIB 
below). The value of was chosen to reproduce the 
viscosity coefficient of a stoichiometric H 2 -air mixture at 
300 K, namely /i//o« 0.2 cm^/s. 

Equations (1) - (8) were solved using a fully unsplit 
corner transport upwind scheme with the PPM spa¬ 
tial reconstruction and the HLLC Riemann solver^®’^^, 
implemented in the code Athena-RFX. This code was 
previously used in a number of turbulent-combustion 
studies^^“^^’^^’^®, which provide further description of the 
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FIG. 1. Overall computational setup in the calculations comprising the NRL Turbulent Premixed Flame Database. Shown is 
the isosurface of the fuel mass fraction, Y, corresponding to the peak reaction rate in case 11 previously discussed in Poludnenko 
et al}‘^ (cf. Fig. 2). 


numerical solver along with the results of detailed con¬ 
vergence tests. 


B. Summary of calculations 

Simulation setup is similar to the one used in our pre¬ 
vious studies of fast turbulent flames.Table II 
summarizes key simulation parameters. 

Calculations are performed in a rectangular domain 
with a uniform Cartesian grid (Fig. 1). The flow is 
initialized with temperature Tq = 293 K and pressure 
Pq = 1.01 X 10® erg/cm®. Kinetic energy is injected at 
the scale of the domain width, L, with a constant rate 
for the duration of a simulation. The amount of kinetic 
energy injected on every time step is ~ 10“® — 10“^ of 
the total kinetic energy.Resulting turbulent flow in the 
upstream cold fuel is homogeneous and isotropic with an 
equilibrium Kolmogorov energy spectrum oc in the 

inertial range extending to scales The method 

of turbulent energy injection used in Athena-RFX (§ II A) 
ensures that in the upstream flow, turbulent integral ve¬ 
locity, C/p and scale, I, are nearly constant both in space 
and time with a standard deviation of < 2% and < 5%, 
respectively. Boundary conditions are periodic in the y 
and 2 directions and zero-order extrapolations in the di¬ 
rection of flame propagation. The absence of unphysical 
effects due to the boundary conditions in our numerical 
approach, in particular, pressure wave reflections from 
the upstream/downstream boundaries, is demonstrated 
below in § IV B. 

Simulation parameters, I and Ui, are chosen primar¬ 
ily to study the dependence of flame dynamics on tur¬ 


bulent intensity (cases S14 - S16 and SlOs) and system 
size (cases S15, S17, S18). In all calculations, I « L/4, 
and the domain width is sufficiently large to ensure that 
L ^ 6l to allow a realistically complex flame structure 
to develop (Fig. 1). 

Combustion regimes surveyed in this work are shown 
in the traditional regime diagram in Fig. 2. Karlovitz 
and Damkohler numbers are®^ 


Ka= — 
tk 




(9) 


7^ '^ed 


( 10 ) 


In the regimes surveyed. Da varies by over an order of 
magnitude (0.1 — 6) and Ka > 1. Note that in this 
work, we define the laminar-flame thickness based on the 
actual flame thermal width, Sl = {Tpy — Tq)/ max(VT), 
rather than the diffusive width, 5 = D{Tq)/Sl- As a 
result, in a stoichiometric H 2 -air mixture (base reaction- 
diffusion model, Table I), Sp^L/^f ~ 52 rather than of 
order unity as is typically assumed.'*’^® Here is the 
kinematic viscosity of cold fuel. This leads to higher Ka 
for a given Ui/Sl compared to the traditional form of 
the regime diagram. 

In all calculations, Sp « 0.032 cm. In all cases with the 
exception of S14 and S18, the internal flame structure is 
resolved with the cell size Aa; = 5^/16, which was pre¬ 
viously shown to provide converged solutions.Two 
lowest intensity cases, S14 and S18, use lower resolu¬ 
tion Aa; = Sl/8. In these cases, the characteristic eddy 
turnover time,^® Ted = l/Ui, is the largest due to lower 
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TABLE 11. Summary of calculations performed 

^ j_ Uj_ IE. Wi WiE 

Ted “ Sl Sl Sl Sl Sl Po Ted Si Si Si 


Calculations that use the base reaction model 


S14 

128^ X 16, 384 

7.07 

3.73 

0.06 

3.73 

1.00 

3.04 

1.65 

5.41 

1.08 

1.25 

0.89 

2.91 

1.32 

S15 

256^ X 16, 384 

22.47 

1.20 

0.06 

3.73 

3.11 

4.90 

2.13 

12.05 

1.21 

2.88 

2.64 

4.36 

2.86 

S16 

256^ X 8,192 

22.46 

0.60 

0.06 

3.73 

6.23 

8.14 

2.09 

21.82 

1.27 

3.28 

5.22 

6.92 

6.00 

S17 

512^ X 16, 384 

9.96 

1.89 

0.06 

7.46 

3.95 

9.99 

3.13 

23.18 

1.48 

1.33 

3.69 

8.27 

4.81 

S18 

512^ X 16, 384 

5.66 

5.99 

0.06 

14.92 

2.49 

9.08 

4.78 

20.75 

1.42 

0.77 

2.35 

8.65 

3.99 





Calculations that use 

a modified reaction model 






SlOs 

256^ X 8,192 

42.09 

0.10 

0.015 

3.73 

38.93 

8.37 

5.16 

16.20 

1.03 

- 

32.66 

33.15 

33.20 

S16s 

256^ X 8,192 

16.20 

0.60 

0.015 

3.73 

6.23 

7.35 

2.10 

14.21 

1.04 

2.81 

5.37 

6.23 

5.58 

S16sa 

256^ X 8, 192 

16.53 

0.60 

0.0075 

3.73 

6.23 

5.56 

2.71 

10.28 

1.01 

2.63 

5.12 

4.84 

5.10 


Calculation that uses 

the base reaction model with temperature-dependent physical viscosity 



S15v 

256^ X 16, 384 

21.36 

1.20 

0.06 

3.73 

3.11 

4.70 

1.65 

11.21 

1.20 

2.67 

2.62 

3.88 

2.62 


See text for the definitions of various quantities. 



FIG. 2. Summary of the combustion regimes modeled in 
the NRL Turbulent Premixed Flame Database. Regimes 10 
and 14 — 18 (green stars) are discussed in this work, while 
other cases were presented in earlier papers, namely case 6 in 
Refs. [12, 13], cases 1 — 11 in Ref. [14], cases 5 — 7 and 12, 13 in 
Refs. [27, 28]. Red diamonds mark cases, in which turbulent 
flames were previously found to exhibit spontaneous DDT.^"* 
Note that DDT did not occur in calculation 10s discussed 
here due to a 4 times lower Sl than in case 10 described in 
Ref. [14]. Solid lines mark Ka = 1 and 100 for Sl = 302 
cm/s, dash-dot lines mark Ka = 1 and 100 for Sl = 75.5 
cm/s (Table I). 


Ui and, in S18, much larger 1. As a result, a com¬ 
pressible hydrodynamic solver, such as the one used in 
Athena-RFX and which is limited by the CFL condition. 


requires a large number of computational time steps to 
advance the solution over one eddy turnover time. This 
significantly increases the computational cost per Ted and 
limits the available resolution. Convergence study per¬ 
formed in case S15, however, showed that at lower [/;, 
the time-averaged turbulent flame speed, St, differs by 
between the resolutions Ax = Sl/8 and 5 l/16. 

Laminar flame speed Sl = 3.02 x 10^ cm/s in calcula¬ 
tions that use the base reaction-diffusion model (Table I), 
i.e., S14 - S18 and S15v. In order to study the effect of 
the laminar-flame Mach number, Mol, values of kq, Dq, 
and B (eqs. 5, 8) in calculations SlOs and S16s were de¬ 
creased by a factor of 4 relative to their values in the base 
model. This resulted in 5^ = 75.5 cm/s (Table I) and 
Mql = 2 X 10“^ vs. 8 X 10“^ in base cases. A lower value 
of S'l in case SlOs also prevented the onset of a sponta¬ 
neous deflagration-to-detonation transition (DDT). Oth¬ 
erwise, it would correspond to case 10 in Ref. [14], which 
had only a relatively short period of quasi-stable flame 
evolution before detonating (Fig. 2). Note that a lower 
Sl implies smaller effective values of ATa for a fixed Ui/Sl 
due to the prefactor {SLE/vfY^'^ (eq. 9). Ka = 1 and 
Ka = 100 for Sl = 302 cm/s and 75.5 cm/s are shown 
with solid and dash-dot lines, respectively, in the regime 
diagram in Fig. 2. 

The effect of the density ratio across the flame, a = 
Po/pp,o, on the flame stability is explored in the calcula¬ 
tion S16sa. Here po and pp^o are, respectively, the fuel 
and product densities. In particular, in this case, the 
value of a was decreased to 3.65 from 7.3 used in all other 
cases, while Sl was kept the same as in case S16s. This 
was done by modifying the values of kq and Dq (eq. 5) 
as well as B, Q, and q (eqs. 3, 8) (see Table I). 

With the exception of case S15v, calculations discussed 
here do not include physical viscosity. Instead, they ex¬ 
plicitly resolve the inertial range of the turbulent energy 
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cascade and use numerical viscosity to provide kinetic- 
energy dissipationd^’^^ It was previously shown^^’^^ that 
this approach accurately captures the flame dynamics in 
high-intensity turbulence, in which the Kolmogorov scale 
rj Sl, or Ka > 1. At the same time, in order to ver¬ 
ify the validity of the results presented here, calculation 
S15v included temperature-dependent physical viscosity 
(eqs. 2, 3, 6 , and 7), with all other simulation param¬ 
eters being the same as in case SI5. In particular, in 
case S15v, Ka « 20, and 77 varied between O.SAx in the 
upstream cold flow and 6.2Aa; in the product, which is 
comparable to the resolution used in other DNS studies 
of high-speed reacting turbulence^°. This calculation is 
discussed in further detail in § IV A. 

Steady-state turbulence is allowed to develop in each 
calculation for « 5Ted- At this point (t = 0), a planar 
unstrained flame is initialized in the y-z plane at the 
distance 6.4L (cases S14, S17, S18) or 7.5L (all other 
cases) from the downstream boundary. After ignition, 
the flame is evolved for another « 7)Ted in order for it to 
become fully developed.Qu anti ties given in Table II, 
including time-averaged values, (...), correspond to times 
t > bTed (or > lOTed since the start of a simulation). The 
total duration of each simulation after t = bred is given 
by ts- In order to achieve total simulation times ts Ted 
(Table II), domains with large aspect ratios ranging from 
32 : 1 to 128 : 1 are used. Simulations are stopped before 
the flame reaches the upstream domain boundary. 


III. RESULTS 

Simulation results summarized in Table II demonstrate 
three key phenomena. First, in all cases studied, St 
is highly variable (also see Fig. 3). Here St is defined 
based on the global fuel consumption rate, m, namely 
St = tn/ pqL'^}^ In case S16, the maximum observed 
value, 5™“^, exceeds St and the minimum observed 
value, S’™*”, respectively, by factors of 2.68 and 10.44. 
Second, in certain regimes, the flame periodically pro¬ 
duces significant overpressures. For instance, in case S17, 
the maximum pressure observed in the domain, 
exceeds Pq by ~ 50%, even though turbulence is highly 
subsonic with Ui corresponding to Mqt = 0.03. Finally, 
the flame can propagate with an average speed, which 
is much larger than the characteristic speed of turbulent 
motions Ui. In case S18, St exceeds Ui by a factor of 
3.65 on average and by a factor of 8.33 instantaneously. 

All three effects are illustrated for base cases S14 - S18 
in Figs. 4 and 5a, which show time histories of normalized 
St/Sl and {P)/Po, along with other quantities discussed 
below. Hereafter, (...) indicates averaging over the flame- 
brush volume, which is bounded by the two y-z planes 
defined as follows. If the x axis points in the direction 
of flame propagation, these two planes have, respectively, 
the maximum (minimum) x coordinate such that all cells 
to the left (right) of the product-side (fuel-side) plane 
have fuel mass fraction Y < 0.05 {Y > 0.95) [12]. 
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FIG. 3. Summary of the turbulent flame speeds in the calcu¬ 
lations discussed in this work. Symbols represent the time- 
averaged turbulent flame speed, St/Sl, in each calculation, 
while the error bars show the corresponding maximum and 
minimum values St^"^/Sl and St'’"'/Sl (see Table II). Note 
that circles correspond to the calculations with the same do¬ 
main width. 


All cases S14 - S18 exhibit pronounced periodic varia¬ 
tions of St, which are also closely followed by the periodic 
increases of (P). While the individual pulses of St are 
irregular, they all have a duration ^Ted- More precisely, 
we define the width of an individual pulse, Tp, as the 
time between two sequential instances when St = St, 
provided that during that time St deviates from St by 
more than 25%. Table H shows that Tp indeed varies 
between 0.77red and 3.28red- Mechanisms of all three 
effects are discussed below. 


IV. MECHANISMS OF FLAME INSTABILITY. 


A. Flame speed variability 


At the values of Ui/Sp considered here, turbulence 
has a minimal effect on the internal structure of a flame 
folded inside the turbulent flame brush.In par¬ 
ticular, average local flame structure reconstructed using 
methods described in Refs. [12, 27] is close to that of a 
laminar flame in all cases. Furthermore, since we consider 
a reacting mixture characterized by Pe = 1 , the local 
burning speed is unaffected by the flame stretch. Conse¬ 
quently, burning occurs in the fiamelet regime, and Sp is 
directly determined by the flame surface area, 


St j At 


( 11 ) 


Hereafter, Ap is defined as the isosurface area of F = 
0.157 corresponding to the peak reaction rate.^^ Ap is 
the global area of the turbulent flame front, or a domain 
cross-section in our case. The factor /c ^ 1 accounts for 
the effect of flame collisions. 
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FIG. 4. Evolution of the normalized turbulent flame speed, St/Sl, (black line, left axis), turbulent velocity fluctuations, 
{SUi}/SL, (green lines, left axis), and cold fuel pressure, {P)/Po, (red line, right axis) in simulations S14 - S17. {SUi} and 
(P) are averages over the flame-brush volume (see text). Horizontal lines show the upstream turbulent integral velocity, Ui, 
(solid), the time-averaged turbulent flame speed, St, (dashed), and the CJ deflagration speed, Scj, (dash-dot) (in cases S14 
and S15, Scj is outside the scale of the graph and is not shown). Note, in panels (a) and (b), the black and red lines are nearly 
coincident. 
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Case S18 
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FIG. 5. Case S18: Evolution of the key dynamical characteristics of the flame, (a) Same as Fig. 4. (b) Normalized, net 
production rates of the vorticity vector components, l^i) jXifi (eq. 16). (c) Normalized individual terms in the vorticity- 

magnitude transport equation, {Itj)/Tifi (eq. 16), for the x and y components of vorticity. Both {li) and {hj) are averages 
over the flame-brush volume (see text). Vertical dashed lines mark the most prominent maxima of St for clarity. Shaded gray 
region in panel (a) corresponds to the time interval shown in Fig. 7b. 


In a turbulent flow, the area of a material surface grows 
exponentially with an e-folding time equal to the charac¬ 
teristic turbulent time scale®’^^’^® 


dAx \ Ax 
dt Jt Ted ' 


( 12 ) 


In a flame, such unbounded growth is counterbalanced 
by two surface-destruction processes. First, flame self¬ 
propagation consumes the surface of a curved front with 
the rate®’^^ 


dAx \ 

dt ) F 


Sl . Sl Ax 


( 13 ) 


where Rc is the characteristic local curvature radius of 
the flame. Second, advection of the flame surface by 


turbulent motions causes individual flame sheets with 
an average separation ~ Rc to collide on a time scale 
~ Rc/Ui [13]. The resulting rate of flame-surface de¬ 
struction is9’38.39 


dAx \ ^ 

dt ) C Rc Rc '^ed 


(14) 


This dynamical picture is frequently modeled in tur¬ 
bulent combustion using various forms of the balance 
equation for the flame surface density.®’^^’^® The relative 
strength of these two flame-surface destruction processes 
is, thus, {dAx/dt)c/{dAx/dt)F = Ui/Sx- 
Based on this, two regimes of turbulence-flame inter¬ 
action can be identified. At low turbulent intensities, 
i.e., at Ui < Slj flame self-propagation is the dominant 
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Case S16s 
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FIG. 6. Case SI6s: Same as Fig. 5. Note that in panel (a), Scj is outside the scale of the graph and is not shown, and black 
and red lines are nearly coincident. Shaded gray region in panel (a) corresponds to the time interval shown in Fig. 7a. 


process of surface consumption. In contrast, in fast tur¬ 
bulence, i.e., when Ui > Sl, destruction of At through 
flame-sheet collisions becomes increasingly more impor¬ 
tant. We emphasize, though, that due to the qualitative 
nature of the above analysis, the boundary between the 
two regimes, i.e., Ui = Sl, is approximate. Furthermore, 
the transition from one regime to another is not abrupt, 
with the relative balance gradually shifting between the 
two flame destruction processes as Ui increases. 

Consider the first, low-speed regime, in which flame 
collisions are rare. Here relatively minimal flame wrin¬ 
kling by turbulence with = DaSL = l{Si/Ui), or 
Rc > I, is sufficient for the time scale associated with the 
flame self-propagation (eq. 13) to become shorter than 
Ted (eq. 12). Self-propagation acts continuously, con¬ 
stantly consuming the flame surface with the rate, which 
varies monotonically with Re- As a result, continuous 
surface creation by turbulence can be counterbalanced by 


the self-propagation. This leads to a quasi-equilibrium 
state, which is characterized by some effective equilib¬ 
rium flame curvature, 

The second, high-speed regime dominated by flame 
collisions differs from the low-speed one in one crucial 
aspect. Unlike self-propagation, flame collisions can be 
highly intermittent, since the flame is consumed only 
when large sections of the flame sheet come into contact. 
As a result, here the balance between surface creation 
(eq. 12) and destruction (eq. 14) can exist only in the 
time-averaged sense. In particular, if at some point Re is 
large (^Oi then {dAT/dt)c <C {dAT/dt)T and At will 
grow exponentially, until eventually surface packing be¬ 
comes sufficiently tight to form an extended region of 
flame collision, which will consume the flame surface. 
Such interplay between the restoring force, which acts 
intermittently, i.e., flame collisions, and the continuously 
acting perturbing force, i.e., turbulence, creates the basis 
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Case SI 6s 
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FIG. 7. Growth of the normalized turbulent-flame surface area, AtIAl, (black lines) and speed, St/SI, (green lines) in the 
course of a typical flame pulsation in simulations S16s (panel a) and S18 (panel b). Red lines show the exponential growth 
of At/Al with the characteristic time scales given in the legend. Time intervals shown in panels (a), (b) correspond to the 
shaded gray regions in Figs. 6a and 5a, respectively. At is the area of the isosurface of the fuel mass fraction Y — 0.157, which 
corresponds to the peak reaction rate. St is normalized by the laminar flame speed, S/, corresponding to the instantaneous 
average fuel temperature and pressure inside the flame brush. 


for the instability. 

This dynamics is illustrated in Fig. 7, which shows the 
evolution of both At/Al and St/Sl (eq. 11) in calcu¬ 
lations S18 and SI6s during the time intervals indicated 
with the shaded gray regions in Figs. 5a and 6a, respec¬ 
tively. After reaching its minimum. At indeed grows 
exponentially as would be expected if only flame-surface 
creation by turbulent stretching were present in accor¬ 
dance with eq. (12). In particular, exponential fits shown 
with red lines are At/Al = 2.01 exp (t/1.55Ted — 10.77) 
in S16s and At/Al = 6.65 exp (t/0.32Ted — 26.88) in S18. 
While in case SI6s the characteristic growth time scale 
is within « 50% of Ted, in case S18 it is a factor of 3 
smaller. The cause of such discrepancy in case S18 will 
be discussed below. 

When At is near its minimum, St/Sl = At/Al al¬ 
most exactly {Ic = 1 in eq. 11). However, as At starts 
to grow, and the flame packing density increases, St/Sl 
begins to deviate from At/Al, i.e., 7e > 1 in eq. 11 
(cf. Fig. 10 in Ref. [20]). Such excess burning rate was 
demonstrated by Poludnenko and Oran to result from 
the formation of cusps caused by flame collisions. In a 
cusp, the volume of the reaction zone is increased by the 
focusing of the heat flux from multiple directions in the 
fuel surrounded by the flame surface. This leads to a 
disproportionately large contribution of cusps to St rel¬ 
ative to the fraction of the overall flame surface area in 
them. The rate of flame-surface consumption in a cusp 
can, in principle, be arbitrarily large, since it is effectively 
the phase velocity, with which two flame sheets collide. 
The burning rate (or fuel consumption rate), however, is 
generally much smaller, though it can typically exceed 
the burning rate of a planar, laminar flame by a factor 


of a few. Both rates depend sensitively on the overall 
geometry of a cusp. When St and At reach their maxi¬ 
mum, the deviation of St/Sl from At/Al is the largest, 
namely 6% in S16s and 22% in S18. At this point, the 
flame collision becomes fully developed. It quickly con¬ 
sumes the flame surface, thus producing the subsequent 
minimum of St and setting up the next pulsation cycle. 
Further detailed discussion of the dynamics of cusps and 
of their effect on the turbulent flame speed can be found 
in Ref. [13]. 

Figure 8 demonstrates the formation of such a flame 
collision in case S16s. Panels (a) and (b) show the flame 
structure at t = 16.72red and t = IQ.lSTe^, respectively, 
i.e., at the beginning and the end of the time interval 
shown in Fig. 7a. During this time, the flame, which 
initially is weakly wrinkled, develops a highly complex, 
convolved structure. Two extended regions are formed, 
in which the flame sheets face each other and are nearly 
in contact (one such region is marked in panel (b), and 
the second one can be seen near the top boundary of the 
domain). Those regions span almost the full width of the 
flame brush. Subsequent rapid flame burn-out in those 
regions restores a much less tangled flame structure, sim¬ 
ilar to the one in panel (a). 

All calculations discussed here have Hi > Sl (Table II). 
Thus, they represent regimes, in which flame collisions 
become progressively more important with increasing Ui, 
according to the analysis above. At the same time, the 
strength of the instability does not increase monotoni- 
cally with Ui, but instead it reaches a maximum at cer¬ 
tain intermediate turbulent intensities (Fig. 3). This is 
demonstrated by cases S14 - S16 as well as SlOs (Fig. 9). 
The variability of St is the smallest in the lowest in- 
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FIG. 8. Case S16s: Change of the flame structure in the course of a typical flame pulsation. Shown is the flame volume bounded 
by the two isosurfaces of the fuel mass fraction Y = 0.05 (red) and Y = 0.95 (blue) at t = 16.72red (panel a) and t = 19.13red 
(panel b). These instances approximately correspond to the times of minimum and maximum flame burning velocity (cf. the 
time period marked with the shaded gray region in Fig. 6a and shown in Fig. 7a). 
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FIG. 9. Case SlOs: Same as Fig. 4. Scj is outside the scale of the graph and is not shown. 


tensity case S14 with the maximum peak-to-peak am- plitude « 3.3. Here, Ui = Sl and flame 
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FIG. 10. Case S16sa: Same as Fig. 4. Note that Scj is outside the scale of the graph and is not shown. Left axis scale is the 
same as in Fig. 6a (the right axis has the same maximum) to facilitate comparison between cases S16s and S16sa. 


self-propagation still has a significant effect (eqs. 13 and 
14), which prevents large growth of At between flame 
collisions. As Ui is increased, variability becomes the 
strongest in case S16, with = 10.4. Fur¬ 

ther increase of Ui/Sl by another factor of 6.2 from case 
S16 to SlOs results only in a marginal increase of St/Sl 
(from 8.14 to 8.37), while decreases effec¬ 

tively back to the value in case S14, namely 3.1. Such 
weakening of the instability in case SlOs is due to a much 
tighter flame packing by the high-speed turbulence with 
the characteristic flame separation i?c —t As a result, 

flame collisions occur more frequently thus interrupting 
periods of the exponential growth of At- In fact, for case 
SlOs, the value of the characteristic pulsation period, Tp, 
is not listed in Table II since pulsations are highly irreg¬ 
ular and their period is not well defined (Fig. 9). Note 
also that in cases S14 - S16, rp/red increases with Ui. 
Such larger rp relative to Ted means longer periods of 
both the exponential growth and subsequent decline of 
At, which, in turn, translates into larger maxima and 
deeper minima of St- 

The effect of the system size, or Z, on the instability is 
more difficult to investigate due to the limited range of 
scales accessible in fully-resolved 3D simulations. Case 
SI7, which has a twice larger I and a somewhat higher 
Ui than case S15 (but the same Ka, or the same tur¬ 
bulent velocity at a fixed scale), shows that 
increases from 5.7 to 7.4 (Fig. 3). In case SI8, which has 
the largest system size, = 4.34. Note, how¬ 

ever, that in this calculation, a somewhat lower Ui rela¬ 
tive to case S15 along with a much larger I translate into 
significantly lower turbulent velocities at a given scale. 

The fluid expansion factor across the flame, a, has a 
strong effect on the magnitude of the flame-speed pulsa¬ 
tions. In particular, in case SI6sa (Fig. 10), a was de¬ 
creased to 3.65 from 7.3 in other calculations. All other 
parameters, including S'p, Sl, and the Zel’dovich number, 
Ze, are the same as in case S16s. Figure 10 shows that 
the overall variability of St is less pronounced in S16sa 
compared with S16s (cf. Fig. 6a), with be¬ 


ing almost twice lower in the former case, namely 3.79 
vs. 6.77 (Table II). 

Finally, as was previously noted in § IIB, all cases 
discussed so far relied on numerical viscosity to pro¬ 
vide kinetic-energy dissipation. A potential effect of the 
temperature-dependent viscosity on the flame dynamics 
described here can be assessed by comparing case S15v, 
which explicitly included physical viscosity (Fig. 11), and 
case S15 (Fig. 4b). Note that both calculations used the 
same base reaction-diffusion model (Table I) and were 
in the same turbulent regime (Table II). The variability 
of St is similar in two cases with « 6.8 in 

S15v compared to 5.7 in S15 and with the average charac¬ 
teristic period of pulsations being almost equal, namely 
Tp = 2.62 and 2.64, respectively. The time-averaged 
St was also very close differing by « 4% (Table II). 
This shows that in high-speed regimes, temperature- 
dependent viscosity does not have a pronounced effect 
on the observed variability of St- 

This result stems from the fact that, as discussed 
above, pulsations of St are driven by flame collisions, 
which are, in turn, formed by turbulent motions on large 
scales > 6l- In high-speed regimes, which are stud¬ 
ied here, scales ^ Sl are part of the turbulent inertial 
range and are, therefore, not sensitive to a specific mech¬ 
anism of kinetic-energy dissipation, i.e., physical or nu¬ 
merical viscosity. For instance, in case S15v, the Taylor 
microscale^^ is A/ = = 0.84(5p, where e is 

the kinetic-energy dissipation rate. At the same time, 
the temperature-dependent viscosity will likely play an 
increasingly larger role at lower Ui as the Kolmogorov 
scale rj ^ Sl and the dissipation range extends to scales 
> Sl, e.g., as in the calculations of Nishiki et al.^^, in 
which r] = 0.65(5l — 0.89(5p. 


B. Pressure variability 

In a previous study^^, we showed that in sufficiently 
high-speed regimes, turbulent flames spontaneously de- 
















13 



FIG. 11. Case S15v: Same as Fig. 4. Note that Scj is outside the scale of the graph and is not shown. All axes have the same 
scale as Fig. 4b to facilitate comparison between cases S15 and S15v. 


velop a runaway process, which results in a rapid pres¬ 
sure increase and may ultimately lead to the formation 
of a detonation. The critical turbulent flame speed was 
shown to equal the speed of a Chapman-Jouguet (CJ) 
deflagration, Scj = Cg/a [3]. In a turbulent flame, Cg is 
the sound speed in cold fuel due to the presence of ex¬ 
tended regions of such cold fuel throughout the interior 
of the flame brush^^ (also cf. Fig. 8b). When St = Scj, 
the amount of energy generated inside the flame-brush 
volume on a characteristic sound-crossing time becomes 
comparable to the thermal energy contained in the same 
volume. As a result, such fast rate of energy release leads 
to the build-up of pressure. 

In an exothermic reaction wave moving with the speed 
S, the ratio S/Scj represents the importance of com¬ 
pressibility effects. This ratio is analogous to the Mach 
number, which characterizes the compressibility of a hy¬ 
drodynamic wave. Hereafter, for brevity, we will refer to 
this ratio as CJ (’’Chapman-Jouguet number”). Note, 
that CJ = aMa. In a typical chemical fuel, a ^ I 
and Scj “C Cg, i.e., CJ deflagrations are highly subsonic. 
In contrast, thermonuclear flames in degenerate plasmas 
have a < 2,^^ and Scj can be close to the speed of sound. 
Values of CJl for a laminar flame in each calculation dis¬ 
cussed here are listed in Table II. 

Figures 4c,d and 5a show that in cases S16 - SIS, St 
(black line) is a significant fraction of Sqj- In particular, 
in S18, on average, CJ = 0.54, while on two occasions 
CJ > I. Pressure (red line) follows St very closely with 
a small lag equal to the sound-crossing time of the flame 
brush, and the overpressures become larger with increase 
in CJ. The fact that P lags behind St shows that it is 
indeed the increase in the burning speed, which drives 
the pressure build-up, rather than vice versa. 

Figure 12 illustrates the evolution of a pressure wave 
in case S18 during one pulsation, the first half of which 
is marked with the shaded gray region in Fig. 5a. In par¬ 
ticular, it shows streamwise, normalized distributions of 
{P)x/Po, where {.■.)x indicates averaging over the domain 
cross-sectional {y — z) planes at each x position. Over¬ 


all, the process of pressure increase inside the turbulent 
flame brush is similar to the one described in Ref. [14] 
(cf. Fig. 3 therein), though peak overpressures observed 
in that work are much higher. A pressure wave begins to 
rise in the interior of the flame brush at x/L ^ 24. As the 
wave grows, it steepens toward the fuel side due to the av¬ 
erage gradient of the sound speed from fuel to product. In 
the spanwise direction, pressure distribution is relatively 
uniform. At t « 8.9Ted, St reaches its peak (cf. Figs. 5a 
and 7b) and starts to decrease due to the flame burn-out 
in the regions of flame collision such as the one shown 
in Fig. 8b. At this point, the pressure wave decouples 
from the flame brush and exits into fuel as a weak shock 
with peak {P)x/Po = 1-42 and Ma = 1.33. A corre¬ 
sponding reverse pressure pulse {{P)x/Po = 1-25) moves 
downstream into product, though it is always weaker and 
broader due to a much higher sound speed in product. 
As both pressure waves (forward and reverse) move away 
from the flame, a rarefaction forms (seen in the last pres¬ 
sure profile at x/L « 24 in Fig. 12), which eventually 
brings pressure and density in the flame brush back close 
to their original values. 

Figure 12 also shows the absence of pressure reflections 
from the open streamwise boundaries of the computa¬ 
tional domain. In particular, dashed black arrows high¬ 
light the propagation of several pressure waves, which 
were formed during previous pulsations. Note that as 
the individual pulses cross the domain boundaries, they 
maintain their shape virtually exactly and no waves mov¬ 
ing in the reverse direction are formed. 

Table II demonstrates that in all calculations consid¬ 
ered, higher average and peak turbulent flame speeds 
result in larger overpressures. In particular, in cases 
S14 - S16, as pulsations become stronger with increasing 
Ui, maximum observed overpressure grows from 8% to 
27% above Pq. System size has a more pronounced ef¬ 
fect. Cases S17 and S18 have larger overpressures (48% 
and 42% above Pq, respectively), even though SJp°‘^/S l 
in them is close to that in case S16. Finally, in cases 
with lower SL^ namely SlOs, S16s, and SI6sa, in which 
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FIG. 12. Case S18: Pressure evolution in the domain in the course of one pulsation. Solid and dashed lines show y — 2 -averaged 
normalized distribntions of pressure, {P)x/Po (left axis), and fuel mass fraction, {Y}x (right axis). Times t = 8.56 — 8.Sired 
correspond to the time period marked with the shaded gray region in Fig. 5a and shown in Fig. 7b. Dashed arrow lines mark 
the propagation of several pressure pulses toward the domain boundaries. Note, for clarity, distributions of {Y)^ are shown 
only for t = 8.56red, 8.94red, and 9.25red. 


CJ = CJl{St/Sl) 1 at all times, overpressures never 
exceed a few percent of Pq^'^ 

For a given system size, a relatively modest increase in 
Ui is required for the flame to transition from a pulsating 
regime, which produces shocks, to a detonation. For in¬ 
stance, case 8 discussed in Ref. [14] (Fig. 2) differed from 
case S16 described here only in terms of Ui, which was 
twice larger. In case 8, after ignition, the flame promptly 
produced a strong pressure runaway and transitioned di¬ 
rectly to a detonation.^® At intermediate intensities, a 
pulsating turbulent flame will exist in a meta-stable state 
characterized by a large probability of developing a pul¬ 
sation of sufficient strength to form a detonation. 


C. Flame self-acceleration 


The third key effect demonstrated by the calculations 
presented here concerns the ability of a turbulent flame to 
propagate with the average speed significantly above the 
characteristic upstream turbulent velocity. This effect is 
most pronounced in the two lowest intensity cases, S14 
and S18, in which St/Ui = 3.04 and 3.65 respectively 
(Table II and Fig. 3). 

Consider the local velocity-fluctuation vector 


SUi{x,y,z) = Ui{x,y,z) - 


T2 


Ui{x,y,z)dydz. (15) 


Here we subtract from the corresponding component of 
the velocity vector its average over the y-z plane. This 


removes the large-scale velocity gradient in the stream- 
wise direction due to the fluid expansion across the tur¬ 
bulent flame. Since we are primarily interested in the 
magnitude of turbulent velocity fluctuations in each di¬ 
rection, hereafter we will consider the absolute value of 
each component of SUi, i.e., SUi = \6Ui\. We then de¬ 
fine the average turbulent velocity fluctuations inside the 
flame brush parallel and transverse to the direction of 
flame propagation, i.e., (SUx) and (SUy^z), respectively. 
Here SUy^z = {SUy + dUz)/2, and (...) again indicates av¬ 
eraging over the flame-brush volume as defined in § HI. 
These velocities allow us to characterize both the average 
intensity and anisotropy of the turbulent flow inside the 
flame brush. Finally, we also define the average upstream 
turbulent velocity fluctuation, SUo, as the average of the 
quantity {SUx + SUy + SUz)/3 over the fuel ahead of the 
flame-brush volume. The definition of SUq combines all 
three components of SUi since the upstream turbulence 
is isotropic. 

Table H shows that in all cases, SUq is close to Ui, 
which is defined based on the kinetic energy spectrum 
in the upstream flow. The difference between the two 
quantities varies in the range « 5 — 20%. In contrast, 
in cases S14, S17, and S18, (SUx) is significantly larger 
than b oth SUp and the average transverse velocity fluctu¬ 
ations, {SUy^z). On the other hand, in these calculations, 
St is much larger than Ui and, at the same ti me, cl ose to 
{SUx)- For instance, in S14 and S18, St and {SUx) differ 
by < 5%. This shows that turbulence inside the flame 
brush undergoes strong anisotropic amplification in the 
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FIG. 13. Case SI8: Conditional pdfs of the turbulent velocity fluctuations transverse and parallel to the direction of flame 
propagation, P{5Uy^z\{Y)x and P{SUx\{Y)z:, (see text for the detailed definition). Pdfs illustrate the anisotropy and inhomo¬ 
geneity of the flow field inside the turbulent flame brush at the times of minimum {t = 8.59red, panels a, b) and maximum 
{t = 8.90red, panels c, d) flame burning velocity in the course of one pulsation (cf. the time period marked with the shaded 
gray region in Fig. 5a and shown in Fig. 7b). Values of {Y}x. in the range 0.01 — 0.99 are split into 16 equal bins. For clarity, 
in each panel, pdfs are shown only for 6 bins, with the legend in panel (a) indicating values of {Y)x. at the center of each bin 
interval. Dashed lines show the pdf of the turbulent velocity fluctuations in the upstream cold flow. Vertical short-dash lines 
mark Ui/Sl, while vertical solid lines mark instantaneous (SUy^z) (panels a, c), and {SUx} (panels b, d). Note, all pdfs are 
normalized by their maximum values to facilitate comparison between pdfs. 


direction of flame propagation, which, in turn, leads to 
the anomalously high flame speeds. 

Complex nature of the change experienced by the tur¬ 
bulent flow inside the flame can be further illustrated 
using conditional probability density functions (pdf) of 
SUx and 6Uy^z, i-e., P{SUj;\{Y)x) and P{SUy^z\{Y)x)- In 
particular, pdfs are conditioned on the average values of 
the fuel mass fraction in the cross-sectional (y-z) planes 
of the domain, (Y)x- Therefore, pdfs, which correspond 
to the decreasing values of (T)x, represent the state of the 
flow field at different positions in the flame brush located 
progressively further from the fuel toward the product. 

For calculation S18, Fig. 13 shows P{5Ux\(Y)x) and 
P{5Uy^z\(Y)x) for two instances, which correspond to the 
beginning and the end of the time interval marked with 
the shaded gray region in Fig. 5a and shown in Fig. 7b. 


Consider first panels (a) and (b), which show the state 
of the turbulent flow field at the beginning of a pulsa¬ 
tion when St is at its minimum {t = 8.59Ted)- At this 
point, PiSUx\{Y)x) and PiSUy^z\{Y)x) at {Y)x = 0.96 
are very close both to each other and to the pdf of the 
velocity fluctuations in the cold upstream flow. In other 
words, near reactants, turbulence inside the flame brush 
is isotropic and is virtually the same as in the upstream 
flow. All three pdfs are Gaussian with the standard de¬ 
viation « Ui (indicated with a vertical dashed line in 
each panel of Fig. 13). As the flow progresses through 
the flame brush toward the products, transverse velocity 
fluctuations undergo only moderate amplification with 
the corresponding pdfs remaining close to Gaussian (cf. 
PiSUy^z\{Y)x) at {Y)x = 0.96 and 0.04 in panel a). In 
contrast, velocity fluctuations in the direction of flame 
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propagation are not only strongly enhanced, but their 
pdfs also become highly irregular (panel b). Note that 
at (Y)x = 0.22 and 0.04, SUx has a nearly uniform prob¬ 
ability distribution over a large range of values. The de¬ 
viation of P{SUx\{Y)x) from the upstream pdf becomes 
particularly large at values of {Y)x 0.5 characteristic 
of the flame reaction zone. 

At the time of maximum St (t = S.Ore^, Fig. 13c,d), 
both SUx and SUy^z near the product side of the flame 
brush at {Y)x = 0.04 become further enhanced compared 
to t = 8.59Ted (red lines in panels c, d). On the other 
hand, at the intermediate values of {Y)x, velocity fluc¬ 
tuations become weaker. For instance, P{SUx\{Y)x) at 
{Y)x = 0.22 is shifted significantly towards lower val¬ 
ues of SUx compared to panel (b), while P{SUx\{Y)x) at 
{Y)x = 0.41 — 0.59 are now close to the pdfs at higher 
fuel mass fractions and in the upstream flow. Transverse 
velocity fluctuations in panel (c) have similar pdfs for all 
values of {Y)x with the exception of 0.04. Note some en¬ 
hancement of both SUy^z and SUx at t = 8.9Ted relative to 
the upstream flow on the reactants side at {Y)x = 0.96. 

Figure 13 shows that in the course of a pulsation, 
turbulence inside the flame brush is not only strongly 
anisotropic, but also pronouncedly inhomogeneous, even 
after accounting for the large-scale velocity gradient. In 
particular, the magnitude of SUx becomes progressively 
larger further downstream, ultimately reaching values, 
which are significantly above the flame-brush average 
(SUx) (vertical solid line in panels b, d). 

The driver responsible for such anisotropic amplihca- 
tion of turbulence can be understood by considering the 
transport equation for the vorticity magnitude, 


' dt 


— LOiUJjSij -j- W, 


Cjjk dp dP 

dxj dxk 


= Xr. 


(16) 

Here, uii is the vorticity vector, Soi = uJilui, Sij is the ve¬ 
locity strain rate tensor, p is density, and Cijk is the Levi- 
Civita symbol. The first, second, and third terms on the 
right-hand side represent, respectively, vorticity produc¬ 
tion by the turbulent velocity strain, effects of compres¬ 
sions or rarefactions, and the baroclinic torque due to 
the misalignment of the pressure and density gradients. 
Since the upstream turbulence is highly subsonic, both 
p and P are nearly constant there, and the dilatational 
and baroclinic terms in the upstream flow are close to 
zero. These terms, however, are important in the flame, 
where p and P can vary substantially. 

The dynamics of all three terms in high-speed turbu¬ 
lent flames was investigated by Hamlington et In 
particular, it was shown that at low turbulent intensities 
when Ui ^ Sl (cf. case FI in Ref. [27] corresponding to 
case 12 in Fig. 2), the dilatational term representing fluid 
expansion across the flame has the dominant effect. In 
contrast, when Ui ^Sl (for a fixed system size, i.e., for 
Da ^ 1), the relative strength of both the dilatational 
and baroclinic terms is minimal compared with the tur¬ 


bulent velocity strain. As a result, in such high-speed 
regimes, the feedback of the flame on the turbulent flow 
field is negligible, and the dynamics of the reacting tur¬ 
bulence approaches that of the non-reacting one. 

In the context of the results discussed here, that analy¬ 
sis had two limitations. First, while the values of Ui in the 
calculations described by Hamlington et spanned a 
similar range as in this work, the domain size, and thus 
I, were twice smaller than the smallest domains consid¬ 
ered here, i.e., in cases S14 - S16 and SlOs (Fig. 2). This 
resulted in lower burning speeds St ^ Scj, and the 
flow did not produce large overpressures observed in the 
present work. Second, the analysis of Hamlington et 
did not consider the potential anisotropy of various terms 
in eq. (16), and instead studied the statistics of the sum 
of each term over all three directions. 

Such anisotropy can arise for two reasons. First, at 
lower C/p the flame-surface normal, and thus Vp, is pref¬ 
erentially oriented in the direction of turbulent flame 
propagation.^^ Second, VP associated with the forma¬ 
tion of pressure pulses is also primarily directed along 
the a;-axis, as discussed above. 

In order to account for such anisotropy, each direc¬ 
tional component of eq. (16) must be analyzed separately. 
The sum of the right-hand-side terms for each direction, 
i.e., for a fixed index i, gives the net production rate, 1^, 
of the magnitude of each vorticity component uji. The 
average value of this rate over the flame brush, {Pi), can 
then be compared with the average of the correspond¬ 
ing net rate in the upstream turbulent flow ahead of the 
flame, 

Figures 5b and 6b show the evolution of {li) for 
each direction in calculations S18 and S16s. In case S16s, 
the average net rate of vorticity production in the flame 
brush is isotropic with values of (Pi) for each direction 
being close at all times. Furthermore, (Pi) are also close 
to the rate in the upstream flow, with the ratio 

of the two varying in the range 0.5 — 1.0. As a result, 
in this case, no excess vorticity is generated inside the 
flame, and there is no net turbulence amplification. This 
is in agreement with the evolution of the average velocity 
fluctuations, {SUi), shown in Fig. 6a. Values of all three 
components of (SUi) are close both to each other and to 
Ui at all times. 

Case S18 exhibits a starkly different behavior. Fig¬ 
ure 5b shows that for the x-component of vorticity, 
{Px)IPx.o is similar to the case S16s and it varies in the 
range 0.4 — 1.4. In contrast, the net production rates 
of the y- and z-components of vorticity inside the flame 
brush are significantly larger than in the upstream tur¬ 
bulence, with {Py)jPyfi and {Pz) jPzfi reaching a value of 
7.1. This results in much higher values of \ujy\ and \uJz\, 
which, in turn, amplifies the velocity fluctuation compo¬ 
nent transverse to uiy and lUz, i-e., SUx- Figure 5a shows 
that {SUx) (solid green line) is indeed The increase 

of {SUx) is partially redistributed among two other veloc¬ 
ity fluctuation components. As a result, {SUy) and {SUz) 
(green dashed and dot-dashed lines in Fig. 5a) are larger 
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than Ui, albeit to a much lesser extent than {5Ux)- 

Comparison of cases S18 and SI6s demonstrates that 
significant increase of the net production rates of the 
transverse vorticity components occurred in the calcu¬ 
lation, which also produced large overpressures, i.e., S18, 
and was absent in the calculation, in which overpressures 
remained at the level of a few percent, i.e., S16s. In 
fact. Fig. 5a,b shows that in case S18 both (ly) and (Iz) 
are closely correlated with the evolution of (P), with the 
peaks of all three quantities being nearly coincident. This 
suggests that increase of {ly) and (Iz) is primarily driven 
by the baroclinic term in eq. (16). 

To demonstrate this. Figs. 5c and 6c show the evo¬ 
lution of the individual terms in eq. (16), ly, averaged 
over the flame-brush volume and normalized by li^. For 
clarity, only terms for the x and y components of the 
vorticity are shown, since terms for the y and z com¬ 
ponents have very similar values. Figure 6c shows that 
in case S16s, the dilatational and baroclinic terms, (I^) 
and (/is), are close to each other and are smaller than 
the turbulent velocity-strain term {In) at all times. Fur¬ 
thermore, terms {Ii 2 ) and (/is) are also similar both in 
the direction of flame propagation and transverse to it. 
Since the dilatational and baroclinic terms terms enter 
eq. (16) with the opposite signs, they effectively cancel 
each other. As a result, the net vorticity production rates 
are primarily determined in this case by the turbulent 
velocity-strain term {In) and the values of {In)/Ii,o and 
{Ii)IIi^o are close. 

In case S18, the balance of all three terms in eq. (16) 
is qualitatively different (Fig. 5c). Here, the dilatational 
and baroclinic terms, {Ii 2 ) and {hs), are significantly 
larger for the transverse {y and z) components of vor¬ 
ticity. Furthermore, the baroclinic term becomes dom¬ 
inant with increasing pressure inside the flame brush, 
as was suggested above. In particular, significant over¬ 
pressures, on one hand, produce large pressure gradients 
(Fig. 12) and, on the other, cause compression and heat¬ 
ing of fuel, which results in a locally thinner flame and, 
thus, in larger density gradients. The combined action 
of these two processes is responsible for the increase in 
the magnitude of the transverse baroclinic term seen in 
Fig. 5c. Note also that in contrast to case S16s, in the 
direction of flame propagation, the dilatational term is 
the dominant one, while the baroclinic term is close to 
the velocity-strain term. 

Close correlation, along with the associated time lag, 
between St, (P), and {6Ux) seen in Fig. 5a demonstrates 
the causal connection between different processes oper¬ 
ating in the course of a flame pulsation. First, growth 
of the flame burning rate (§ IV A) leads after a short 
time lag to the build-up of pressure (§ IV B), formation 
of large pressure gradients, and increase of the relative 
strength of the baroclinic term in eq. (16). Subsequently, 
after a somewhat longer delay, resulting higher produc¬ 
tion rates of the transverse vorticity components amplify 
the turbulent velocity in the direction of flame propaga¬ 
tion. This causal connection is further corroborated by 


case S16s (Fig. 6). It shows that the exponential growth 
of St does not require either significant pressure build-up 
or anisotropic turbulence amplification inside the flame 
brush. Therefore growth of St is the root cause, rather 
than the consequence, of the latter two processes. 

In general, larger overpressures do not directly trans¬ 
late into a larger increase of the turbulence intensity in¬ 
side the flame. Instead, the strength of this turbulence 
generation process depends sensitively on the relative 
balance between the turbulent strain and the baroclinic 
terms in eq. (16). Consider cases S14 and S16. Even 
though peak pressures produced in S16 are significantly 
larger than in S14, St, {5Ux), and {SUy^z) are all close 
to Ui- I n cont rast, in S14, an 8% peak overpressure re¬ 
sults in {5Ux) a factor of ^3 larger than Ui. This shows 
that with increasing Ui, the relative strength of the tur¬ 
bulent strain term in eq. (16) grows rapidly allowing it 
to dominate the baroclinic torque even in the presence 
of 30% overpressures. In case SlOs, in the limit of very 
high turbulent intensities and in the absence of pressure 
pulsations, turbulent velocity fluctuations are virtually 
identical both inside the flame brush and in the upstream 
flow (Table 11). 

Finally, this process of turbulence amplification allows 
us to explain the characteristic time scale of the expo¬ 
nential growth of the flame surface area in case S18. 
In particular, it was discussed above that in Fig. 7b 
this time scale is a factor of three shorter than Ted, 
namely Q.?)2Ted- Table II, however, shows that in this 
case {SUx)/Ui = 3.47. Therefore, the eddy turnover time 
scale corresponding to {SUx) is = 0.29Ted, which is 
within 10% of the time scale in Fig. 7b. This shows that 
high turbulent intensity inside the flame brush acceler¬ 
ates the exponential growth of At and decreases the char¬ 
acteristic growth time scale. Furthermore, such faster 
turbulence also results in a shorter average period of pul¬ 
sations, Tp, as can be seen in cases S14, S17, and S18 
compared to other calculations (Table II). 


V. DISCUSSION AND CONCLUSIONS 

We presented results of a suite of fully resolved cal¬ 
culations modeling the interaction of a steadily driven 
turbulence with a premixed flame. These calculations 
span a broad range of regimes characterized by Ui > Sl 
and Da = 0.1 — 6.0. 

One of the key conclusions of this work is that in 
certain regimes turbulent premixed flames are intrinsi¬ 
cally unstable. In particular, even under the most ideal¬ 
ized conditions of a statistically steady, homogeneous, 
isotropic, Kolmogorov-type upstream turbulence, the 
flame burning speed can exhibit periodic pulsations with 
the observed peak-to-peak amplitude > 10. 

Such values are significantly above the peak-to-peak am¬ 
plitudes of < 2 — 3, which were previously found in the 
calculations of turbulent premixed flames done by us^^’^^ 
as well as by other groups^^’^^. Note that this is, in ef- 
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feet, a limit-cycle instability in contrast with the flame 
instability described by Poludnenko et which re¬ 

sults in the unbounded growth of St and, ultimately, in 
the transition to a detonation. 

In the course of pulsations, the flame speed can become 
a large fraction of, or even exceed, the speed of a CJ de¬ 
flagration, Scj- This results in the formation of pressure 
pulses. With increasing turbulent intensity, such pulses 
become progressively stronger forming shocks, and, ul¬ 
timately, they can produce a detonation. The figure of 
merit characterizing the strength of such compressibility 
effects in a flame is the ratio St/Scj = CJ (’’Chapman- 
Jouguet number”), rather than the Mach number. We 
emphasize that this process does not require any con- 
hnement of the flow, i.e., presence of walls or obstacles, 
and it can occur in a highly subsonic upstream turbu¬ 
lence. In particular, formation of shocks was observed 
with Ma > 1.3 in the course of flame interaction with 
the Mar = 0.02 upstream turbulence. 

Results presented here also show importance of the 
baroclinic torque in fast turbulent flames. Overpressures 
produced in the course of pulsations can couple with the 
density gradients across the flame through the baroclinic 
term in the vorticity evolution equation and result in an 
efficient mechanism of vorticity generation. Even rela¬ 
tively small overpressures of < 10% (case S14) can am¬ 
plify turbulent intensity inside the flame brush signifi¬ 
cantly above the levels in the upstream flow and increase 
the overall turbulent burning velocity. This process is, 
effectively, a reacting-flow analog of the mechanism un¬ 
derlying the Richtmyer-Meshkov instability.^® The differ¬ 
ence is that the pressure gradient is not imposed exter¬ 
nally, e.g., through a shock impacting a perturbed inter¬ 
face with a density jump. Instead, pressure is generated 
by the perturbed interface itself, i.e., by the turbulent 
flame. 

It is important to emphasize the difference between 
the process of turbulence amplification discussed here 
and the classical picture of flame-generated turbulence, 
which traces its origins back to the early experimental 
work by Karlovitz et In both cases, the underlying 
physical mechanism of turbulence generation is the same, 
i.e., coupling of pressure and density gradients inside the 
flame through the baroclinic torque^®’^®’"^®. The origin 
of pressure gradients, however, is qualitatively different. 
Traditionally, flame-generated turbulence is considered 
in the context of lower-speed flames, which experience a 
pressure drop associated with the fluid expansion across 
a flame. In contrast, fast turbulent flames discussed here, 
which propagate at near- or super-CJ speeds, produce a 
pressure rise. In realistic chemical flames, the magnitude 
of a pressure drop, and thus of the associated pressure 
gradients, is typically relatively small compared to the 
overpressures, which can be produced in super-CJ defla¬ 
grations. Such overpressures can, in principle, reach a 
significant fraction of the von Neumann pressure in a CJ 
detonation before the resulting shock waves would cause 
a spontaneous DDT. Consequently, high-speed turbulent 
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FIG. 14. Normalized average turbulent flame speeds, St/Sl, 
in each calculation shown as a function of the normalized av¬ 
erage turbule nt velo city fluctuations in the direction of flame 
propagation, {SUx)/Sl (see Table II). Note, both axes have 
the same scale as Fig. 3 to facilitate comparison. 


flames can produce much higher levels of flame-generated 
turbulence, which, in turn, would have a significantly 
larger effect on the overall burning velocity. 

Results discussed here suggest several directions, in 
which further development of existing turbulent com¬ 
bustion models is warranted.First, large variabil¬ 
ity of the turbulent flame speed implies that in high¬ 
speed regimes, the state of a turbulent flame cannot be 
characterized by a single value of Instead, a 

given upstream turbulent intensity, Ui, and integral scale, 
Z, must be associated with a probability distribution of 
flame speeds, which can span a large range of values. 

Second, St is primarily controlled by the turbulent ve¬ 
locity fluctuations in the direction of flame propagation, 
6Ux, inside the flame brush rather than in the upstream 
flow. In fact, the large scatter of values of St seen in 
Fig. 3 is removed, when St is plotted as a function of 
(SUx) rather than [//. In particular, Fig. 14 shows that 
all cases, with the exception of the highest-intensity case 
SlOs, cluster tightly around the line 

^=Sl+W^- (17) 

This shows that a flame can exhibit strong self¬ 
acceleration, when 6Ux is significantly amplified inside 
the flame brush. Therefore, it is important for turbulent 
combustion models to predict accurately the magnitude 
of such anisotropic turbulence amplification Eventually, 
however, with the increase in turbulent intensity, eq. (17) 
breaks down, as evidenced by case SlOs (see § IV A; 
Ref. [13] provides further discussion of the flame dynam¬ 
ics in such extremely high-speed regimes). 

Third, in the regimes, in which the flame is suscepti¬ 
ble to large-amplitude pulsations, St depends not only 
on the turbulent-flame surface area, At, but also on the 
details of the flame configuration. In particular, peri¬ 
odic formation of extended regions of flame collisions re¬ 
sults in the excess burning velocity, which is represented 
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by the proportionality factor Ic in eq. (11) (cf. Fig. 7; 
also see Ref. [13]). This is the case even in the absence 
of any flame-stretch effects as in Le = 1 reacting mix¬ 
tures. Our ability to predict accurately the magnitude of 
Ic requires deeper understanding of the properties of such 
flame collisions. Over the years, several studies have con¬ 
sidered their effect on the flame dynamics.How¬ 
ever, a number of basic questions remains unanswered, 
e.g., the rate of formation of flame cusps on different 
scales, their morphology, effects of turbulent intensity 
and system size, to name a few. 

Finally, another important consequence of turbulence 
amplification inside the flame brush for turbulent com¬ 
bustion modeling is that this process can cause the flame 
to evolve in a combustion regime, which is different from 
the one determined based on the upstream turbulent con¬ 
ditions, e.g., thin reaction zones rather than the corru¬ 
gated flamelets regime.®^ In fact, in certain cases as a 
result of flame self-acceleration, we observed DDT, even 
though upstream turbulent intensities were nominally 
not high enough for the detonation to form. These cal¬ 
culations will be discussed in a separate paper. 

The fundamental underlying mechanism, which drives 
the variability of the flame speed, is the interplay between 
turbulence, which continuously acts to create the flame 
surface, and the intermittently occurring flame collisions, 
which provide the restoring mechanism consuming the 
flame surface. While this picture provides a qualitative 
explanation of the flame dynamics that was observed in 
the calculations presented here, a number of important 
questions remains. 

The key question concerns our ability to predict the 
amplitude of pulsations. In particular, it is possible to 
predict very accurately the rate of growth of the flame 
surface (cf. Fig. 7). At the same time, it is not clear 
what controls the length of time, during which At under¬ 
goes the exponential growth. Figure 7a shows that such 
growth can last several large-scale eddy turnover times 
required to form flame collisions on the largest scales in 
the flow. This suggests that some process suppresses for¬ 
mation of flame collisions on smaller scales, which evolve 
on much shorter time scales. 

One possible explanation for this suppression is the 
redistribution of turbulent kinetic energy in the flame 
between different scales. It was discussed by Hamling- 
ton et that fluid expansion across the flame results 
in an anisotropic suppression of vorticity, which weakens 
turbulent motions perpendicular to the flame surface on 
small-scales and enhances them on large scales. Effec¬ 
tively, eddies, which rotate in the plane perpendicular to 
the flame surface, expand as they pass through the flame. 
It is, however, those eddies that fold the flame and pro¬ 
duce flame collisions. Suppression of small-scale eddies 
by fluid expansion would prevent the formation of flame 
collisions on small scales and would allow the exponential 
growth of At over long time periods. 

Such explanation agrees with the observed weakening 
of the instability when the fluid expansion factor, a, is 


decreased. This situation is analogous to the Landau- 
Darrieus instability^®, for which a (along with the Mark- 
stein number) determines the smallest unstable wave¬ 
length. It is not clear, however, whether in the case of 
the flame instability discussed here, large-amplitude pul¬ 
sations become impossible below a certain a, or if the 
decrease in a can be offset by a larger system size. 

Elucidating the effect of fluid expansion is important in 
order to understand whether thermonuclear flames in de¬ 
generate plasmas are also intrinsically unstable. On one 
hand, thermonuclear flames typically have a < 2 [43]. 
On the other, the turbulent integral scale during the lat¬ 
ter stages of a thermonuclear Type la supernova explo¬ 
sion is orders of magnitude larger than the laminar flame 
width.®® Large variability of the flame burning speed, 
as well as the resulting formation of pressure pulses and 
flame self-acceleration, would have a significant effect on 
the dynamics of Type la supernovae compared to the 
current explosion models®’^^’®® 

A related question concerns the overall dependence of 
the instability on the system size. In particular, it is 
not clear whether the amplitude of flame-speed pulsa¬ 
tions can become larger with increasing 1. It was dis¬ 
cussed above that flame collisions are important in the 
regimes, in which Ui > Sl- However, this condition can 
be satisfied in large systems, 5l, simultaneously with 
Da ^ 1. Such large values of Da imply that flame self¬ 
propagation cannot be neglected. Thus, these situations 
may in fact represent a hybrid regime, in which flame self¬ 
propagation is dominant on smaller scales, while flame 
collisions form on large scales allowing the instability to 
operate. 

It is important to understand the connection of the 
flame instability discussed here with the traditional ther¬ 
moacoustic combustion instabilities.^’®’^ In particular, it 
is not clear how the major mode of the pulsating in¬ 
stability, which has a period ~ Ted, would couple with 
the acoustic field of a combustor, which can have a dif¬ 
ferent mode imposed by the combustor geometry. Fur¬ 
thermore, it was recently suggested that laminar flames 
may be intrinsically unstable even in the anechoic com¬ 
bustion chambers, i.e., in the absence of any resonant 
interaction with the acoustic field inside a combustor (T. 
Poinsot (2014), private communication). Similarly to the 
dynamics described here, such flames also exhibit peri¬ 
ods of significant flame-surface growth as they are being 
stretched by the flow, which are followed by rapid flame 
burn-out when the flame forms a highly cusped structure. 
It would be important to understand a possible connec¬ 
tion between that intrinsic instability and the instability 
described here. 

Finally, effects discussed in this work can have impor¬ 
tant implications in practical combustion applications. 
Both the variability of the burning speed and flame 
self-acceleration through the baroclinic torque can place 
a combustion system outside of its designed operating 
regime. Furthermore, the ability of a subsonic turbulent 
flame to generate strong pressure pulses and shock waves 
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can have a destructive effect on a combustion system. In 
this context, an interesting question concerns the con¬ 
nection of the instability described here with the engine 
knock, which is typically associated with the propaga¬ 
tion of spontaneous reaction waves.On the other hand, 
controlled propagation of flames in this unstable regime 
can provide a distinct mechanism of pressure-gain com¬ 
bustion. This regime would, in effect, be intermediate 
between the isobaric turbulent combustion, e.g., used in 
the traditional turbine engines, and the detonation-based 
combustion used in pulsed and rotating detonation en¬ 
gines. 
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